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Abstract 

The lattice instabilities of perovskite structure oxides are responsible for 
many of their interesting properties, such as temperature-dependent ferro- 
electric phase transitions. First-principles calculations using linear response 
theory provide an accurate means to determine the lattice dynamics through- 
out the entire Brillouin zone (BZ). Using the LAPW linear response method 
Q, we have previously carried out such a study on ferroelectric KNb03 [||. 
We present here the results of a similar investigation for cubic SrTi03. While 
KNbC>3 has only ferroelectric-type instabilities, SrTiOa exhibits both ferro- 
electric (FE) and antiferrodistortive (AFD) instabilities. We correctly predict 
the known instability at the R-point in the Brillouin zone, which is responsi- 
ble for the AFD phase transition to the tetragonal structure at about 105 K. 
Furthermore, the phase space of the ferroelectric instability is greatly reduced 
compared to KNb03. Antiferrodistortive instabilities exist in one-dimensional 
cylindrical tubes extending along the entire R-M-R line in the Brillouin zone. 
The essentially 1-D character of these tubes corresponds to real-space planar 

instabilities characterized by rotations of oxygen octahedra. 
PACS numbers: 63.10.+a, 77.80.-e, 77.84.Dy 
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I. INTRODUCTION 



Strontium titanate has been the subject of ongoing theoretical and experimental investi- 
gation for over 30 years, due to its unusual properties. In cubic perovskite structure KNb0 3 
and BaTiC>3, transitions to ferroelectric phases occur as the temperature is lowered. SrTiOs, 
which has the same high-temperature structure, behaves somewhat differently in that it un- 
dergoes an antiferrodistortive (AFD) transition as the temperature drops below about 105 K 
P|. This is due to a softening of the lowest frequency triply degenerate zone-corner R-point 
vibrational mode of lowest energy MM. ^ n this mode, the Sr and Ti atoms remain fixed, 
while the oxygen octahedra rotate about one of the cubic axes passing through the Ti atom. 
The sense of the rotation is opposite in adjacent cells in all directions as shown in Fig. [l|. 
As the temperature continues to drop below the AFD transition, there is a sharp rise in 
the value of the static dielectric constant. This increase appears to follow a Curie- Weiss 
law consistent with an impending ferroelectric (FE) transition at about 20 K. With decreas- 
ing temperature, however, the static dielectric constant saturates at a high value, and the 
material does not become ferroelectric even at the lowest temperatures |6|]7[]. It has been 
suggested that the FE transition is inhibited due to quantum paraelectric behavior, in which 
the zero-point motion of the atoms suppresses the long-range FE order ||||. It is thus of 
great interest to characterize the structural instabilities of SrTiOs. The present work is con- 
cerned with the wavevector dependence of structural instabilities in SrTiOs, obtained from 
first-principles calculations, and their interpretation in light of previous similar calculations 
on KNb0 3 . 



II. METHODS 



Our calculations are carried out utilizing the LAPW formalism ||10|| , wherein the 
wave functions, charge density, and potential all have dual representations. Within non- 
overlapping muffin-tin spheres, these quantities are expanded as numerical radial functions 
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times spherical harmonics, while in the interstitual region, a plane wave representation is 
used. This method has several attractive features, such as employing an unbiased representa- 
tion in the interstitual region and the ability to easily treat first-row and d-band atoms. The 
major drawback presented by this method is the use of position-dependent basis functions. 
This introduces correction terms when computing such quantities as the forces on the atoms 
|rj|l , and thus complicates the computer programs. In order to investigate all possible soft 



phonon modes in a particular substance, we must map out the dispersion curves throughout 
the Brillouin zone. This is accomplished quite easily using linear response theory |Ti[T2|. In 
linear response theory, the total variational energy is expanded to second order as a function 
of the equilibrium charge density and the first-order change in the charge density: 

E T = [p (0) ] + \EM [p (0) ] + X 2 E^ + ■■■. (1) 

Minimizing this energy with respect to the first-order change in the charge density yields a 
pair of equations which must be solved self-consistently: 

N 

p (1) = E^ (1) , (2) 

1=1 

(^ o) -4 o) )# ) =-(k ( ;}[p (i) ])^ (o) . (3) 

The equations are the first-order analogs of the usual local density approximation (LDA) 
equations, i.e. this is an exact theory; there are no approximations beyond LDA. 

The dynamical matrix elements for a given wavevector q are found via the first-order 
forces. The forces on all the atoms are calculated for at most 3N selected atomic displace- 
ments at each wavevector, where N is the number of atoms in the primitive cell. To map the 
phonon dispersion curves throughout the Brillouin zone, we compute the dynamical matrix 
on a uniform grid of q-points. The short-range real-space force constants are then obtained 
by Fourier transform, after subtracting out the analytic long-range dipolar contribution. By 
inverse Fourier transform, adding in the long-range part, we may compute the dynamical 
matrix at an arbitrary q-point. For the SrTiC>3 calculations, we have used a 6x6x6 mesh of 



q-points. Exploiting symmetry, we had to compute the dynamical matrix at only 20 q-points 
in the irreducible BZ wedge. Typically, at each q-point, 3N=15 self-consistent calculations 
were performed. 

The electronic contribution to the static dielectric constant, e^, is determined from the 
induced potential: 

V ind (q) 



q • £oo • q = lim 



(4) 



y total (q 

The macroscopic polarization at long wavelength is given by the first order change in the 
charge density 



P = lim - p (1) (q), 
q^o q 

and from this we can compute the Born effective charges 



(5) 



zLpU) = Zj + adP a /dT P {j) 



E=0 



(6) 



Here, Z, is the ionic contribution, and the second term is the unit cell volume times the 



J i 3 



derivative of the polarization with respect to displacements of the atoms at zero electric 
field. 



III. RESULTS 



We first performed self-consistent total-energy calculations in the cubic phase to deter- 
mine the equilibrium properties. Table | lists the theoretical lattice parameters and bulk 



moduli for the present work, along with those of King-Smith and Vanderbilt ||13|| , as well 
as some experimental values. Our calculated lattice parameter is only 0.4% larger than the 
experimental value, so we have used the theoretical value in the following calculations. The 
electronic component of the static dielectric constant, ^oo, was found to be 6.63, which is 
about 28% larger than an experimental value of 5.18. This kind of overestimation of is 
typical of LDA calculations P^JTSfl . 
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Table |TT| lists the Born effective charges (Z*) for each type of atom. For Sr and Ti, the 
values are isotropic, while for the O atoms, there are two distinctive values — for displace- 
ments along and perpendicular to the Ti-0 bonds, labeled || and _L, respectively. Also listed 
are the results of calculations by Zhong, King-Smith, and Vanderbilt Respective theo- 
retical lattice parameters were used for the calculations. Note the anomalously large values 
for Z*(Ti) and Z*(0)y due to strong covalent interactions between these atoms |17|,18| . We 



find that our results are generally in good agreement with those of Zhong et al. ||16|| . 

Table |T| compares our zone-center optic phonon frequencies with planewave pseudopo- 
tential frozen phonon calculations and experimental data. Both calculations find unstable 
transverse optic (TO) modes at the T-point with imaginary frequencies. The longitudinal 
optic (LO) mode frequencies were obtained using 

d lo d to | 47re 2 (Zf ■ q) a (ZJ ■ q), 

where D TO is the zone-center dynamical matrix without macroscopic field, Z* and Mj are 
the Born effective charge tensor and mass for atom i, Q is the volume of the unit cell, a, (3 
are Cartesian indices, and q is a unit wavevector. The results of Zhong, King-Smith, and 
Vanderbilt were determined using £00=5.18 from experiment, whereas we used our larger 
calculated value of 6.63. To determine what effect £00 has on the LO mode eigenvalues, a 
second calculation was performed using £00=5.18. The highest LO mode is most sensitive, 
changing from 751 to 832 cm~l, while the other two are fairly insensitive. 

We present the calculated phonon dispersion curves for cubic SrTi03 plotted along high- 
symmetry directions in Fig. |2|. The TX, TM, and TR lines correspond to the (100), (110), and 
(111) directions, respectively. Imaginary frequencies are represented by negative values. The 
character of the modes at the zone-center and zone-boundaries has been labeled according 
to the notation by Cowley |19j . The form of some of our lower dispersion curves for the TM 
and TR directions is roughly in accordance with data reported by Shirane and Yamada at 
120 K ||, including an apparent softening of a mode at the R-point. Similar behavior is 
seen at the M-point, but the lowest zone-boundary phonon mode M3 was reported as being 



temperature-independent through the AFD transition. 

Since we are mainly concerned with structural instabilities and how they relate to the 
observed phase transitions, we will focus on the portions of the dispersion curves which lie 
below the cc;=0 dashed line in Fig. |2|. Note the large phase space for unstable modes. Our 
calculations indicate unstable modes at the R-point (R25), zone-center (r 15 ), and M-point 
(M 3 ). These instabilities are of two types: FE (r^), and AFD (R25 and M 3 ). In the Ti 5 TO 
mode, Ti atoms move parallel to one of the Ti-0 bonds, and the oxygen octahedra move 
in the opposite direction. This mode is responsible for the FE transitions in materials like 
BaTiOs and KNbOs [0,[n[] • The M3 mode is nearly identical to the R25 mode in Fig. [l] with 



one exception. The rotation of the octahedra is in the same sense in neighboring cells along 
the vertical c-axis, but remains opposite in the horizontally adjacent cells. The dispersion 
curves of Fig. ^| show |co>(r 15 )| < |a;(i2 2 5)|, suggesting that LDA may find a FE structure 
(rhombohedral?) to be lower in energy than the AFD tetragonal structure. Of course 
experimentally there is no observed transition from AFD to FE phase, although it is seen 
in the Monte Carlo simulations of Zhong and Vanderbilt [ 22] . Recently, they have included 



quantum fluctuation effects into these simulations, and find the FE transition completely 
suppressed . 



The regions of instability in the BZ are better visualized in Fig. ||], in which isosurfaces, 
corresponding to uj=0, are shown for the lowest unstable phonon modes. The cubic BZ, 
with T at the center, is also shown. The inner isosurface is centered at the r-point, and 
can be visualized as three interpenetrating "cookies", one perpendicular to each Cartesian 
direction. The region of FE-type instability is interior to this first isosurface. Between the 
first and the second isosurface, which lies near the zone edges, all modes are stable. Thus 
unstable modes are present along the entire R-M-R edge of the BZ. In a repeated-zone 
scheme, this isosurface would appear as narrow cylindrical tubes extending from R-M-R. 

We now discuss the current work in light of results for similar calculations on cubic 
KNbOs by Yu and Krakauer || . The dispersion curves of KNbOs also reveal a large region 
of instability for the lowest modes. In contrast to SrTiOs, a soft mode is present at the 
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X-point, and the lowest R-point mode is stable. Although M-point instabilities exist for 
both materials, the character of these modes is quite different. In KNb03 the motion of the 
atoms is the same as for the zone-center Tis mode with the exception that two O atoms 
perpedicular to the motion remain fixed, i.e. it is essentially the "ferroelectric" soft mode. 
In SrTiOa, the motion is a rotation of the octahedra as described above. A similar uj=0 
isosurface plot for KNbC>3 shows that the FE instability extends from the zone-center all 
the way to the zone-boundary within planar slab-like regions. An AFD instablility was 
not present in the KNb03 calculations. Notice that in SrTiC>3, the large planar regions of 
FE instability have shrunk down to three interpenetrating "cookies", and the phase space 
corresponding to the FE character is greatly reduced. 

Because the regions of instability in KNbOs reside in slabs which extend across the 
entire BZ, any linear combination of phonon modes within one of these regions will also be 
unstable. In a real-space picture, unstable atomic chain-like structures can develop, where 
the Nb atoms are all displaced in the same direction The width of the slab-like regions 

in KNbOs indicates the length of the shortest such chains to be about 5 lattice parameters 
0. We can perform a similar qualitative analysis of the real-space motion of the atoms in 
SrTiOa from its isosurface plot. Since there exists a continuous instability extending from R- 
M-R, linear combinations of phonon modes along this cylindrical tube-like region will also be 
unstable. The character of modes along the zone edges is one of rotating oxygen octahedra; 
thus, an unstable thin disk or planar region could form in real space wherein the octahedra 
are rotated in opposite senses in neighboring cells, with the surrounding regions remaining 
undisplaced. From the width of the cylinders containing the instability, we estimate the 
smallest radius of the real-space disks to be approximately 3-5 lattice parameters. Because 
of the sensitive dependence of the soft modes on volume, it is possible that these cylinders 
may pinch off at the M-point at slightly reduced volumes. We have not investigated this 
yet. 
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IV. SUMMARY 



In summary, we performed first-principles linear response calculations of the wavevector 
dependence of structural instabilities for cubic SrTi0 3 . In comparison with previous work on 
KNbC>3, we notice three differences. First, calculations on KNbC>3 show a ferroelectric- type 
instability only, while SrTi0 3 calculations exhibit both FE and AFD instabilities. Second, 
the phase space of the two-dimensional FE instability (corresponding to chains in real space) 
is greatly reduced and no longer extends to the zone-edges as in KNbC>3, but is restricted 
to the interior of a region likened to three inter-penetrating "cookies". Third, the AFD 
instability occupies a rather large region of phase space extending along the entire BZ edges 
in cylindrical tubes. The essentially 1-D character of these tubes corresponds to real-space 
planar instabilities characterized by rotations of the oxygen octahedra. Future calculations 
for SrTi0 3 in the tetragonal phase (below the AFD transition) will be required in order to 
investigate any instabilities — especially ferroelectric instabilities — that persist. 
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FIGURES 

FIG. 1. Motion of the oxygen atoms in the zone-corner R25 phonon mode. The Ti atoms lie 
inside the octahedra and do not move. The Sr atoms are also fixed but are not shown. The 
octahedra rotate about the c-axis in opposite senses in all adjacent cells. 



FIG. 2. Calculated phonon dispersion curves for cubic SrTiOs at the theoretical lattice constant. 

FIG. 3. Zero-frequency isosurfaces of the lowest unstable phonon modes over the entire BZ. 
The T-point is located at the center of the cube. Unstable modes exist inside the central surface 
and along the full length of the zone-edges. 
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TABLES 



TABLE I. Equilibrium properties determined by calculations and experiment for SrTiOa 



Lattice parameter a 


Bulk modulus 


(a.u.) 


(GPa) 



Present 
PW 
Exp. fc 
Exp. c 



7.412 
7.303 
7.380 



190 
200 

179±6 



Tlanewave ultrasoft-pseudopotenial calculation by King-Smith and Vanderbilt fl3 
6 Wyckoff, Crystal Structures, vol. 2, 2nd ed. (1964). 
C T. Mitsui et ai, Landolt-Bornstein Series, Group III, vol. 3 (1969). 



TABLE II. Born effective charges for cubic SrTi03. 





Present 


PW-BP a 


Lattice parameter (a.u.) 


7.412 


7.303 


Z*(Sr) 


2.55 


2.54 


Z*(Ti) 


7.56 


7.12 


Z*(0)|| 


-5.92 


-5.66 


Z*(0) ± 


-2.12 


-2.00 



a Planewave Berry phase calculation by Zhong et al. [16 . 
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TABLE III. Calculated zone-center optic phonon frequencies (in cm 1 ) in cubic SrTi03. 



Present 



PW a 



Exp. 6 (300 K) 



TO modes: 

lOOi 

151 

219 

522 

LO modes: 
146 A 146 s 
439 A 449 s 
751 A 832 B 



41t 
165 

546 

158 
454 
829 



Soft 
175 

IR-silent 
545 

171 

474 
795 



"Reference |16[ . 
^Reference fE$. 

^ Using Eq.(7) and the calculated value £oo=6.63. 



B 



Using Eq.(7) and £oo=5.18 extracted from experiment. 
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This figure "fig3.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/9703240vl 



